Second Quantized Reduced Bloch 
Equations and the Exact Solutions for 
Pairing Hamiltonian 

Liqiang Wei 

Institute for Theoretical Atomic, Molecular and Optical Physics 
Harvard University, Cambridge, MA 02138 

Chiachung Sun 
Institute of Theoretical Chemistry, Jilin University 
Changchun, Jilin 130023 P. R. China 

February 2, 2008 

Abstract 

In this article, we present a set of hierarchy Bloch equations for 
the reduced density operators in either canonical or grand canonical 
ensembles in the occupation number representation. They provide a 
convenient tool for studying the equilibrium quantum statistical me- 
chanics for some model systems. As an example of their applications, 
we solve the equations for the model system with a pairing Hamilto- 
nian. With the aid of its symplectic group symmetry, we obtain the 
statistical reduced density matrices with different orders. As a special 
instance for the solutions, we also get the reduced density matrices of 
the ground state for a superconductor. 
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1 Introduction 



The language of the second quantization is a convenient and ubiquitous tool 
for studying the quantum many-body systems in modern physics. The quan- 
tum statistics of the identical many particles is automatically satisfied by 
the commutation relations of the creation or annihilation operators, and the 
complicated mathematical manipulation with wavefunctions such as Slater 
determinants and with matrix elements of operators is significantly simpli- 
fied. It is no surprise that a familiarity with the related knowledge of the 
occupation number representation is a must for a serious theorist. 

In a recent paper, we have developed an equation of motion approach 
for a direct determination of the reduced density matrices for the statistical 
ensembles under the framework of equilibrium quantum statistical mechanics 
[1-4]. We have derived the hierarchy Bloch equations which provide the laws 
according to which the reduced statistical density operators vary in temper- 
ature. We also take the orbital approximation, and solve the equations for 
the identical fermion systems with two-body interaction in a grand canon- 
ical ensemble. We have gained significant new physical insights. We find 
that the usual Fermi-Dirac statistics for the free fermion gas holds also for 
the interacting fermions, and that the usual Hartree-Fock equation at zero- 
temperature can be extended to the case of any finite temperature. The 
mean force field that every orbital or particle is subject to will be from the 
contribution of both coherent and incoherent superpositions of the wavefunc- 
tions of other orbitals or particles, and is therefore temperature-dependent 
[5-8]. This opens a way for investigating the interplay between the micro- 
scopic structure and macroscopic thermodynamic quantities for molecules 
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and solids [5-9]. In addition, the Fermi-Dirac distribution for interacting 
particles demonstrates that there exists an effective one-particle description 
for the interacting many-body system [2,10]. 

In the current paper, we begin the presentation of the scheme for an exact 
solution of the reduced Bloch equations which will be developed in a second 
quantized form. In quantum many-body theory, the systems are often repre- 
sented by the model Hamiltonians in an occupation number representation. 
We study the pairing Hamiltonian, which is a useful model for the study of 
nuclear structure [11,12], microscopic state of superconductors [13,14], and 
electronic structure of molecules [15,16]. 

We organize this paper as follows. In the next Section, starting from the 
Bloch equation, we deduce a set of hierarchy equations in the occupation 
number representation for the reduced density operators in the canonical 
ensemble. The derivation is based on the second quantized form of the con- 
traction operator which is presented in the Appendix. We also extend it 
to the case of a grand canonical ensemble. In Section 3, we first briefly re- 
view the symplectic group symmetry of the pairing Hamiltonian and present 
the second quantized form of its eigenvectors characterized by the particle 
number, the weight and the seniority number. Then using the symmetric 
conditions among D 1 , D 2 and D 3 , we solve the reduced Bloch equations in 
the canonical ensemble for the pairing Hamiltonian and obtain its canonical 
reduced density matrices with different orders. As a special situation, they 
yield the reduced density matrices of the ground state for the superconductor. 
We conclude our discussion in the final Section. 
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2 Reduced Bloch Equations in the Occupa- 
tion Number Representation 

We consider an identical iV-fermion system with a two-body interaction. Its 
Hamiltonian has the following form, 

# = X>(i) + X>(i,j), (i) 

i=l i<j 

where h(i) and v(i,j) are one- and two-particle operators, respectively. As- 
sume a generic and complete set of single-particle states {\i >}, which form 
a single-particle Hilbert space V. For each \i >, the corresponding creation 
or annihilation operators are defined as 

\i >= af\0 >, <n\j >= 6ij\0 > . (2) 

They satisfy the anticommutation relations, 

{a;,a+} = 5ij, 
{ai,aj} = {a+,a+} = 0. (3) 

Define an iV-fold antisymmetrized direct product space of the single-particle 
Hilbert space V as 

V NA — ViAV 2 A ... A V N . (4) 
The Fock space F can be written as 

oo 

F=J2 ®V NA , (5) 

N=0 

and any operator T on F can be expressed as the direct sum of the iV-fold 
tensor T N , 

oo 

T=J2 ®T N . (6) 

N=0 



The contraction operator, which defines the reduced density matrix as follows 
[1,17,18] 

D p — L P N {D N ), (7) 

is a superoperator acting on the iVth-order density matrix D N , and has 
following second quantized form (see the Appendix), 

L P N {T N ) = ^(A-) N - p (T N ), (8) 

where A_ is defined by 

A-(T N )=Y,< TN W- (9) 

i 

Equipped with these definitions as well as the Hamiltonian in the occu- 
pation number representation, 

H = Y1 hijdfaj + Vijkiaj alj aia k , (10) 

ij ijkl 

where hij =< i\h\j > and v^i =< ij\v\kl >, we can now deduce the second 
quantized reduced Bloch equations [19-25]. For the canonical ensemble, the 
Bloch equation is [1-4] 

- = ill hijafaj + i Vijkiata+cna^j D N , (11) 

where the D N is regarded as an iVth-order irreducible tensor, and (3 is the 
reverse of the product of the Boltzmann constant ks and absolute temper- 
ature T. Applying on both sides of above equation with the operator L P N 
gives 

- jL DP = E hijL p N (ala,D N ) +W v ijkl L% (afafa^D") . (12) 

° ' ij ~ ijkl 



The evaluation of the first term of the right hand side of the Eq. (12) results 
in 

^L P N (aja ]D N ) = (N — p)j^^ ( D N ) af + ^aja 3 L p N ( D N ) , 

(13) 

and the calculation of the second term yields 

m rv( + + n N\ (N- P )(N-p-l)N\ j-p+2/nJVx + + , 

—L p N [aJaJaia k D J = (p + 2)! a i a k L N (D jafaj + 

+ {N (p7l)\ l { a ! a ^ L N +1 (D N )af - aia iak If N + \D* 
+ ^aja1-aia k L p N (D N ). 

Therefore, the set of reduced Bloch equations for the canonical ensemble in 
the occupation number representation is 

^ " \ ij ijkl 

Define 

= 7;(%fc; - Vjiki), (16) 
then the above equation can be simplified as 

- ^-D p = HD P + ^ZZ E h l]aj D p+l at + E %«a+a iflfc ^ +1 a+ + 

For the grand canonical ensemble, whose pth-order reduced density ma- 
trix is defined by [1] 

D P G =t®( N v )L P N [D N ], (18) 

N=p \ P / 



similar procedure leads to the following equation 




hd g + h ij a i D G +la t + Vijkidj 'a x a k D V Q l af + 



ij ijkl 



+ 9 2 VijkiaidkD^ 2 af aj , 



(19) 



ijkl 



where 




(20) 



and 



H = hijafaj + - ^ Vi jk iafajaia k . 



(21) 



ij ijkl 



This is a set of second quantized reduced Bloch equations for the reduced 
density matrices in the grand canonical ensemble. 

3 Reduced Density Matrices with Different 
Orders for Pairing Hamiltonian 

3.1 The Pairing Hamiltonian and its Eigenvalues and 
Eigenstates 

Prior to solve the reduced Bloch equations in the case of the canonical en- 
semble for a pairing Hamiltonian, we need to review its eigenvalue problem 
as well as related symmetry properties which is discussed in details in paper 



[26]. 



Label a set of 2A spin-orbitals > (dim V = 2 A = r) as 



1, i, A, 1 , % , A 
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where % — i + A, and \i > represents the pairing of a single-particle state 
\i>. The pairing Hamiltonian is defined as [26-29] 

H = eN-GQ+Q- 

r A 

= tJ2 a t a i- G Yl a t a s a j' a J> ( 22 ) 

i=l i,j=l 

where 

i=i 

and 

A 

Q- = &i = a *' a *' (23) 

i=i 

are called the quasi-spin operators. They are the generators of the Lie algebra 
for the quasi-spin group SU®(2). The second term of the Hamiltonian indi- 
cates that the interaction is only for the pairing states of the particles, and 
is therefore called the pairing force. The matrix element of the interaction 
is represented by a constant G. When G > 0, the interaction is repulsive; 
when G < 0, the interaction is attractive. 

Since the pairing Hamiltonian Eq. (22) exhibits a symplectic symme- 
try, its iV-particle eigenstates should be the basis vectors of the irreducible 
representation space of the symplectic group, and can be characterized by 
its weight [30]. For the symplectic group Sp(2X) and its corresponding Lie 
algebra C\, the weight is the eigenvalue for the generator, 

Su = a+ai- ata^, i=l,2,...,A, 

which is the difference of occupation numbers between the pairing spin or- 
bitals \i > and \i >, with the value being 1, 0, or —1. They can be grouped 
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together as a vector 

(S) = (si,s 2 ,-,sa). 

Furthermore, every irreducible representation can be characterized by the 
maximal weight. For the symplectic group Sp(2X), this maximal weight takes 
the form of (l u X ~ U ), where the number v is called the seniority number, 
and the irreducible representation of the symplectic group can be written 
as < V >. However, there is a possibility that the same weight might 
still correspond to the different basis vectors. To completely distinguish all 
the basis vectors of an irreducible representation space for the symplectic 
group, we therefore need to classify the iV-particle eigenstates according to 
the following group chain 

Sp{2\) D Sp{2\ - 2) D ... D Sp{2). (24) 

This means that each basis vector will further be characterized by a group 
of A seniority numbers: 

(V) = (z/i l/ 2 ... V\). 

Thus, the eigenvectors of the pairing Hamiltonian will be completely iden- 
tified by the particle number N, the weight (S), and seniority number (V): 
\N, (S), (V) >. They belong to the irreducible representation spaces < V 1 >, 
< V 2 >,..., and < V x > of the group Sp(2), Sp(4), and Sp(2X), respec- 
tively, and satisfy the eigenequation 

S U \N,(S),(V) >= Sl \N,(S),(V) >, ( Si = 1,0,-1). (25) 

The branching rule of V AA under the group chain (24) of S p (2\) is 

[1 N ] ^<1 N > + < 1 N - 2 > +...+ < 1 > or < >, 
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for N < X, or 
[1*] 



>< l 2 ^ > + < l 2X - N - 2 > +...+ < 1 > or < > (26) 



for N > A. The dimension of the irreducible space < V > is given by 

dim <l»>=( 2X )-( , 2 \ ) ■ (27) 



The analytical expression for the basis vectors can be constructed as the 
following form [26] 



\N,(S),(V) >= \N,< V >,W >-- 



(X — v — m)\ ^ m , , 

(28) 



where 
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{k-v k + 2)(k-v k + l) 



k-l 



j2bt-(k-u k + i)bt 



,i=i 



(29) 



and 



m=-(N-u). 



In Eq. (28), the ki is the positive number which gives the difference of two 
consecutive seniority numbers in the (V) which is two: v k — v k -i — 2; the ji 
means the same thing but its value equals one: v\j\ — = 1. The sign of 

ji is decided by the sign of sy\: if sy\ — 1, j > 0; while if sy\ = —1, j < 0. 
In other words, ±k±, ...,±k y are the pairing numbers in the corresponding 
symplectic table W, while j±, ...,j x are the unpairing numbers appearing in 
the table, which captures all the features of the basis vectors (28) classified 
according to the group chain (24) [26,31,32]. 

The actions of the quasi-spin operators Q + or Q_ on the eigenvectors are 



Q + \N, (S), (V) >= -y (N — v + 2)(r — N — v)\N + 2, (S), (V) >, (30) 
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or 



g_|7V, (S), (V) >= y( N -u)(r-N-u + 2)\N - 2, 0), (V) >, (31) 

respectively. Therefore, we can get the analytic expression for the energy 
eigenvalue as follows 

E(u) = Ne- ^(N -v)(r- N -is + 2). (32) 



3.2 The Exact Solutions of Reduced Bloch Equation 

We are now in a position to solve the reduced Bloch equation (15) for the 
pairing Hamiltonian (22). For the case of p = 1, the set of hierarchy equations 
can be written as 

-^D 1 = NeD 1 -\{N -l)G^Q-D 2 at + \{N - l)G£a+Q_I> 2 aJ 
(7V-l)(7V-2)^ n3 



GQ-D A Q+. (33) 
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Since the pairing Hamiltonian is invariant under the symplectic group Sp(2X), 
the density matrix describing the states of the system can be characterized 
by its irreducible representations as follows 

D N (P) = J2e- mu) D N (u), (34) 

where 

D N {v)= I Yl \N,<r>,W><N,<l v >,W\ (35) 

dim < V > w ^ lv> 

and 

TrD N {y) = \ {v = N, N - 2, 1 or 0). (36) 
11 



In addition, 

*Ee-e E M = Z(J3,v,N). (37) 

For the same reason, the pth-order reduced density matrix of iVth-order 
density matrix can be written as 

D*>(P) = Y,u; p (v,v,(3)Di>(v), (38) 

where 

D p (u')= J2 \p,<l"' >,W><p,<l v ' >,W\ (39) 
w&<i v ' > 

{y =p,p-2,...,l or 0) 

and 

TrD p ({3) = Z(P,V,N). (40) 

The thermal probability coefficients u> p (is, u ,/3) are to be decided. 

According to the branching rule (26) and dimension formula (27), the one- 
particle Hilbert space V (denoted as < 1 >) is irreducible and its dimension 
is r. Thus 

D 1 (P)=J2M"A,{3)D 1 (1)- (41) 

V 

The two-particle Hilbert Space V 2A comprises two irreducible subspaces < 
> and < l 2 > of Sp(2X). Their dimensions are 1 and ^ ^ ^ — 1> respectively. 
Therefore 

D\(3) = £ 0, (3)D 2 (0) + £ tJ2(v, 2, (3)D 2 (2). (42) 

V V 

The three-particle Hilbert space V 3A also includes two irreducible subspaces 
< 1 > and < l 3 > of Sp(2X) with dimensions being r and 
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respectively. Thus 

D\(3) = J2 u 3 (v, 1, (3)D\l) + J2 3, f3)D 3 (3). (43) 

V V 

The substitution of the expressions (41), (42), and (43) into Eq. (33) then 
leads to an equation that the above five coefficients satisfy. Together with 
three trace conditions (40) of D 1 , D 2 and D 3 as well as the iV-representability 
conditions among D 1 , D 2 and D 3 as shown below 

D 2 {(3) = Lj[D 3 ((3) 



J2"3(vA,P)Ll [D 3 (l)\ +J2^3,f3)Ll [D 3 (3)\ , (44) 



and 



D 3 (3) , (45) 



D\?) = Ll{D 3 ((3) 

= 5>3(M,/?)4[£> 3 (1)] +E^3(^,3,/3)4 

v v 

it seems that we can completely determine the reduced density matrices 
for the pairing Hamiltonian. However, due to the complexity of the iV- 
represent ability, it is often replaced by other conditions. Here, we choose the 
symmetry conditions among D 1 , D 2 and D 3 . 
Suppose 



L\ \D 3 (l) =a 1 D 2 (0) + b 1 D 2 (2). 



From 



Q+Q- = A 



N 



r»(\gxg\), 



(46) 



(47) 



where is the expansion operator (see references [17,31] or the Appendix), 
we get 



ai 



Tr 



\g X g\Ll(D 3 (l)) 



3<r-2). 



(48) 
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But 



TrL\ 



D\l)\ =oi + 6i-(r + l)(r-2) = r, 



we therefore have 



61 = 



3(r-2)' 



In the same manner, we assume 



(49) 



(50) 



T 2 



D 3 (3) = a 2 D 2 (0) + b 2 D 2 (2), 



and obtain 



a 2 = 0, b 2 



r(r — 4) 



3(r-2) 

The insertion of Eqs. (46) and (51) into Eq. (44) yields 



(51) 



(52) 



D 2 ((3) = 5>3( M ,/?) ! — ^ 2 (0)+£ 



^3(1/, 1, P)^-^ + W3(i/, 3, /9)^-4ll D 2 (2). 



3(r-2) 



3(r- 2) 
(53) 



Similarly, 



Lj(D s (l)) = ^(1), 

4(^(3)) = i(r + l)(r- 4)^(1), 



and we get the expression for -D 1 (/3) 



1 



^M, /?) + 3, %(r + l)(r - 4) 

6 



(54) 



(55) 



Therefore, with the aid of symmetry conditions, we have expressed the ther- 
mal probability coefficients of D l (f3) and D 2 (f3) by those of D 3 ((3). 

Now we can calculate three terms of the right hand side of Eq. (33) as 
follows, 



^ d i=i 



(56) 
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Y,atQ-D 2 (P)at = - ]>> 3 (^, 1,0)-^- £ \i >< i\, (57) 

i 



r - 2 r 

v Z i=l 



Assume 



w 3 (v,l,/3) 



-PE{v) 



W 3 (f, 1), 



(59) 



where 



exp(-/3E(y)), v = N, N - 2, 1, 
are linearly independent. Then from the trace condition of D 3 (f3), 



u 3 (u,l)r + u 3 (u,3)-(r + l)(r-A) 



= £ 



we get 



uj 3 (u, l)r + w 3 (f, 3)-(r + l)(r - 4) = I. 



(60) 



(61) 



Furthermore, the substitution of Eqs. (56), (57) and (58) into Eq. (33) gives 



£■ 



.^- >)(f - ff - y + a)+tt(M)G ^-'Mr-') 



= 0. 
(62) 



This leads to 



- j(N - v){r - N - v + 2) + uj 3 (u, 1)G ^- ^ 2) - 0. (63) 
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Hence 



and 



co 3 (u, 1) = 



3(N-v)(r- N-u + 2) 
N(N-l)r(r-2) 



lo 3 (u,3) 



6 

r(r — 4)(r + 1) 



3(iV-z/)(r-iV-z/ + 2) 
N(N-l)(r-2) 



(64) 
(65) 
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Finally, we obtain the first-order, the second-order, and the third-order re- 
duced density matrices of iV-particle systems governed by the pairing Hamil- 
tonian at any finite temperature as shown below, 



D\(3) 



l j2\i><i\Z(P,V,N), 



r 7=i 



(66) 



^ = Ee-^ ( ^-^ 1) 7 + 2) ^(0) + 



X 



1 - 



_ mv) f A{N-v){t-N-v + 2) 

\ N(N - l)r(r - 2) 2 (r-2)(r + l) 

3(iV-z/)(r-iV-/y + 2) 
N(N-l)(r-2) 



D 2 (2), 



(67) 



and 

£ 3 (/?) 



V e -^M 3(iV ~ ^ (r ~ - ~ - + 2) I) + V e -^M ^ 



X 



N(N — l)r(r — 2) 
3(7V-z/)(r-7V-i/ + 2) 



D 3 (3). 



r(r — 4)(r + 1) 
(68) 



7V(7V- l)(r-2) 

These density matrices decide the microscopic states of the system, from 
which we can calculate the thermodynamic properties. 



4 Conclusions 

We have derived the second quantized form of the reduced Bloch equations 
for both the canonical and grand canonical ensembles. The deduction is 
based on the occupation number representation of the contraction operator, 
and is similar to the ones we used before as described in papers [1] and [2]. 
The equations obtained provide a direct route for exactly solving the reduced 
density matrices for the systems represented by the model Hamiltonian in 
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the occupation number representation. They thereby open an avenue for 
investigating the microscopic structures of these systems and their relations 
to the equilibrium thermodynamic properties. 

As an example for exactly solving the reduced Bloch equations, we have 
obtained the reduced density matrices of the pairing Hamiltonian with dif- 
ferent orders as shown in Eqs. (66), (67) and (68). One very interesting 
conclusion we can make is that these reduced density matrices can yield the 
projection of the ground state of a superconductor in the iV-particle space 
V N . This is the case when G < and N is an even number, and the ground 
state of the paring Hamiltonian (22) is an AEGP (Antisymmetrized Extreme 
Geminal Power) function | N, (0 A ), (0 A ) >. Therefore, if we choose E(0) as 
the energy zero point and let u — and [3 — > oo in Eqs. (66), (67) and 
(68), we can obtain the reduced density matrices of the ground state for the 
superconductor as follows 

D\T = 0) = ^El<X<| = ^i, 

D^T-0) - r - N + 2 DH0) + 2{N - 2) DH2) 

r-N 2r(N-2) . . 

\g >< g\ + t \,, AT '— D 1 A D , 



(r-2)(N-l)^ ai (r-2)(N-iy 



and 



D ^-^ - (N-l)r{r-2) D Uj r(r - 2)(r - 4)(7V - 1)^ [6) 

9(r-N) . 



r(r - 2)(r - A)(N - 1) 



(r-2)(r-4)(JV-l) v ; 



These results are identical to Coleman's [33]. 
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Since there are many other interesting physical systems often modelled 
in an approximate way by the Hamiltonians with a second quantized form, 
we expect that the ideas and methodologies we have discussed in this paper 
can be further explored, and that the fundamental equations (15) or (19) 
we have derived will find wide applications in the study of their microscopic 
states as well as their interplay with macroscopic thermodynamic properties. 



A Appendix: Second Quantized Form of Con- 
traction and Expansion Operators 



Consider an antisymmetrized iV-particle state ket |iii2...ijv > composed of 
the single-particle states defined in the Section 2. It is easy to see that the 
action of an annihilation operator cij on this state ket is 



where the superscript " indicates that the state \ik > is destroyed. The con- 
jugate of Eq. (70) is 

N 

< i 1 i 2 .-.i N \aj~ = ^2(-l) N ~ k 5 jik < iii 2 -.-i k ...i N \. (71) 
k=i 

The combination of Eqs. (70) and (71) therefore yields 



A? 



aj\ixi 2 ...i N >= ^2(-l) N k 6 jik \i 1 i 2 ...i k ...i N >, 
k=i 



(70) 



E(-i) fc+fc s t 



> ikj ,\lll2--lk-lN >< 3lJl-3 k '--jN\- 



(72) 



18 



On the other hand, the state ket |iii 2 ...ijv > can be expanded as a deter- 
minant as follows, 



1 N 



|iii 2 ...ijv >= -7= XX -1 )* 1 |»i*2--»fc---*jv > |«* >(at) • (73) 
viV fc=1 

Thus, the action of contraction operator on the basis operator \iii 2 ...iN >< 
3i32--3n\ is 

■N-l 



1 

TV 



1 ' 

= TFE(- 1 ) fc+fc< W i i i 2-"**-»*Jv >< hh-k'-M- ( 74 ) 



1 (|zii 2 ---«Af >< Jij2---jAr|) = -T7j2 a i\ i ^-" i N >< 3i32---3n\(4 ■ (75) 



The comparison of Eqs. (71) and (74) leads to 

1 

N 

Furthermore, since any iVth-order tensor operator T N can be expressed as a 
linear combination of the basis operators {\i\i2---iN >< 3i32---3n\}-> we finally 
get 

L%-\T N ) = ^ ai T N 4, (76) 

i 

which is the second quantized form for the contraction operator. 

The expansion operator is the reverse operation of the contraction 
operator L P N and transforms a joth-order tensor to an iVth-order one (p < N) 
[17] 

T^(T p ) = A N [T p <g> I N ~ P ] A N , (77) 

where is the antisymmetrized operator for iV particles. Since the ex- 
pansion operator is the adjoint mapping of the contraction operators L P N 
[17] 

< r^(T p ), T N >=< T p , L P N {T N ) >, (78) 
19 



we can also obtain its second quantized form as follows 

^iM^E^-V (79) 

i 

References 

[1] L. Wei and C. C. Sun, Physica A, in press (2003); e-print, 
|cond-m at/0306306, 

[2] L. Wei and C. C. Sun, Physica A, in press (2003); e-print, 
|cond-m at/0306307| 

[3] F. Bloch, Zeits. f. Physick 74, 295 (1932). 

[4] J. G. Kirkwood, Phys. Rev. 44, 31 (1933). 

[5] K. Pichler, D. A. Halliday, D. D. C. Bradley, P. L. Burn, R. H. Friend, 
and A. B. Holmes, J. Phys.: Condens. Matter 5, 7155 (1993). 

[6] W. Klemperer, Ann. Rev. Phys. Chem. 46, 1 (1995). 

[7] J. Yu, M. Hayashi, S. H. Lin, K.-K. Liang, J. H. Hsu, W. S. Fann, C.-I. 
Chao, K.-R. Chuang, S.-A. Chen, Synth. Met. 82, 159 (1996). 

[8] S. P. Kwasniewski, J. P. Francois, and M. S. Deleuze, J. Phys. Chem. A 
107, 5168 (2003). 



20 



[9] H. G. Drickamer and C. W. Franck, Electronic Transitions and the High 
Pressure Chemistry and Physics of Solids (London, Chapman and Hall, 
New York, 1973). 



[10; 

[11 
[12 

[is; 

[14 
[15 

[ie; 

[17; 
[is; 

[19 

[20; 

[21 
[22 
[23 



W. Kohn and L. J. Sham, Phys. Rev. A 140, 1133 (1965). 

(a) B. Lorzao, Nucl. Phys. A397, 225 (1983); (b) B. Lorzao and C. 
Quesne, Nucl. Phys. A403, 327 (1983). 

K. Hagino and G. F. Bertsch, Nucl. Phys. A679, 163 (2000). 

J. M. Blatt, Prog. Theor. Phys. 23, 447 (1960). 

R. Balian, H. Flocard and M. Veneroni, Phys. Rep. 317, 252 (1999). 

B. Weiner and O. Goscinski, Phys. Rev. A 27, 57 (1983). 

B. Weiner, H. J. Aa. Jensen and Y. Ohrn, J. Chem. Phys. 80, 2009 
(1984). 

H. Kummer, J. Math. Phys. 8, 2063 (1967). 

J. E. Harriman, Phys. Rev. A 17, 1257 (1978). 

L. Cohen and C. Frishberg, Phys. Rev. A 13, 927 (1976). 

H. Nakatsuji, Phys. Rev. A 14, 41 (1976). 

H. Schlosser, Phys. Rev. A 15, 1349 (1977). 

J. E. Harriman, Phys. Rev. A 19, 1893 (1979). 

C. Valdemoro, Phys. Rev. A 45, 4462 (1992). 

21 



[24] D. A. Mazziotti, Phys. Rev. A 57, 4219 (1998). 

[25] A. J. Coleman, Intern. J. Quantum Chem. 85, 196 (2001). 

[26] Z. H. Zeng, C. C. Sun, and A. J. Coleman, in Density Matrices and 
Density Functionals, pp. 141-165 (1987), eds by R. Erdahl and V. H. 
Smith Jr. (Dordrecht: Reidel). 

[27] S. G. Nilsson and O. Prior, Mat. Fys. Medd. Dan. Vid. Selsk. 32, no 16 
(1960). 

[28] R. W. Richardson, J. Math. Phys. 18, 1802 (1977). 

[29] F. Pan and J. P. Draayer, Phys. Lett. B 442, 7 (1998). 

[30] B. G. Wybourne, Classical Groups for Physicists (Wiley, New York, 
1974). 

[31] Z. H. Zeng and C. C. Sun, J. Mol. Sci. 4, 287 (1986). 
[32] A. J. Coleman, J. Math. Phys. 27, 1933 (1986). 
[33] A. J. Coleman, J. Math. Phys. 6, 1425 (1965). 



22 



